Strong Coupling and Double Gap Density of States in Superconducting MgB 2 
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Using scanning tunneling spectroscopy at T = 4.2 K, we perform simultaneously the topographic 
imaging and the quasiparticle density of states (DOS) mapping in granular MgB2. We observe a 
new type of spectrum, showing a pronounced double gap, with the magnitudes of Ag = 3.9 meV and 
Al = 7.5 meV, i.e. well below and well above the BCS limit. The largest gap value gives the ratio 
2A_t/fcsT c = 4.5, which implies strong electron-phonon coupling. Other superconducting regions 
are found to have a characteristic BCS-shaped DOS. However, the variation of the spectral shape 
and lower gap widths, from 2.0 meV to 6.5 meV, indicate the importance of surface inhomogeneity 
and proximity effects in previously published tunneling data. Our finding gives no evidence for any 
important gap anisotropy. Instead, it strongly supports the multiple gap scenario in MgB2 in the 
clean limit, and the single gap scenario in the dirty limit. 

PACS numbers: 74.20.-z, 74.50.+r, 74.70.Ad 



Since its announcement by Nagamatsu et al. in Jan- 
uary 2001 p], the discovery of superconductivity in the 
simple metallic compound MgB2 has stimulated a great 
effort in the field. A very high critical temperature 
Tc — 39 K, well above the values reported for all other 
conventional superconductors, and quite simple binary 
structure, make this material a good candidate for tech- 
nological applications. At the same time this unusually 
high Tq opens a new fascinating question about the elec- 
tronic structure and pairing mechanism involved. 

Numerous experimental reports confirm that MgB2 is 
a BCS-type phonon-mediated superconductor. A sig- 
nificant boron isotope effect ( a ~ 0.26 ) was observed 
in both magnetization and specific heat measurements 
pj. Measurements on nuclear spin-lattice relaxion rate 
in the n B NMR study by Kotegawa et al. || have shown 
the existence of an s-wave superconducting state with a 
large isotropic gap (2A/ksTc — 5) well above the weak- 
coupling BCS value of 3.52. The strong coupling scenario 
could be explained considering the high phonon frequen- 
cies (uj ~ 700 K) in MgB2 due to the light boron mass 
M. To complicate the picture, it has been shown theoret- 
ically that MgB2 has a particular band structure in which 
both two- and three-dimensional bands are present and 
that, in the clean limit, could even give rise to a multiple 
gap superconductivity ||. 

The direct measurement of the superconducting DOS 
of MgB2 is urgently needed. Photoemission experiments 
have already been performed but are, unfortunately, 
limited by energy resolution. Details of the magnitude 
and symmetry of the superconducting gap can be ob- 
tained with tunneling spectroscopy (TS), giving directly 
the quasiparticle DOS with an energy resolution of a few 
kT. If the full scanning TS (STS) is performed, the spa- 
tial mapping of the local DOS is possible. It has also 
proved useful in the case of inhomogeneous systems 
Finally, in the tunneling conductance, spectral features 



found above the gap allow a self-consistent analysis of 
the electron-phonon interaction B . 

Unfortunately, the published tunneling spectra clearly 
diverge |I(Jjl4|] and no STS data have yet been presented. 
While the shapes of the spectra indicate an isotropic gap 
for the pairing symmetry, the gap widths are quite dif- 
ferent, all being too small to be realistic for the bulk 
MgB2. Moreover, in jll| the spectra were found ex- 
ceedingly smeared. To explain these different results, 
even a very sophisticated symmetry was suggested for an 
anisotropic order parameter [[PI . In this series tunneling 
experiments |Io| - |l4|| no evidence whatever for a multiple 
(or double) gap scenario was shown. Consequently, an 
important uncertainty exists in establishing the shape of 
the quasiparticle DOS and the gap magnitude pertinent 
to the bulk superconductivity. 

In this Letter we show that the complete STS, even on 
the small scale of a few hundred nanometers, allows one 
to interpret the wide variation in the previous tunneling 
data. The measurements were performed on MgB2 pow- 
der (Tc = 38.7 K), which was glued to the sample holder 
by silver paint and mechanically etched Pt/Ir wires were 
used as STM tips. Measured locally, the conductance 
spectra reveal a distinct double gap structure as shown in 
Fig.l. These show little or no smearing other than ther- 
mal broadening at 4.2 K. As we shall argue, this double 
gap is an intrinsic property of MgB2 and not a proximity 
effect @. 

The tunneling conductance in Fig.l is also character- 
ized by a very flat spectral background, good sharpness, 
conservation of quasiparticle states, and by no shape de- 
pendence on the tip-sample distance. These conditions 
demonstrate the good quality of the tunneling junction. 
The spectra show a relatively small gap ~ 5 mV with al- 
most no excitations inside. It is followed by two further 
peaks, that are seen in both occupied and empty states 
sides. As we shall describe below, these features corre- 
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spond to the larger second superconducting gap A = 7.5 
meV. It exceeds all those reported in the tunneling lit- 
erature ]lO|-[l^|. Consequently, the corresponding ratio 
2A/fcsTc = 4.5 clearly points to strong electron-phonon 
coupling. 
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FIG.l : Double gap structure in local conductance spectra 
(dots). A small gap As = 3.9 meV and a large one Al = 7.5 
meV exist together. Flat spectral background and zero states 
at the Fermi energy are characteristic features of such spec- 
tra. Solid line: A generalized BCS fit, using Eqs. (1-3), with 
the parameters F = meV and spectral weights Cs = 0.94 
and C L = 1 - C s = 0.06. 

An analogous double-gap spectrum was very recently 
observed in an inverted SIN junction JlTj and in an inde- 
pendent point-contact spectroscopy experiment p^ |. ] 
has also been shown that both gaps close at the sam 
temperature |L7|, i.e. the critical temperature Tc of tb 
bulk material. Thus both gaps are intimately related t 
the superconducting properties of MgB2. This behavic 
is consistent with the theoretical prediction of Liu et a 
H for two gap superconductivity in the clean limit fc 
bulk MgB 2 . 

A double gap, such as in Fig.l, is remarkable but coul 
a priori originate from different physical scenarios. I 
the first one, the small gap reflects a proximity supei 
conductivity induced in a thin metallic surface layer. I 
this case, the observed bump indicates the position c 
the superconducting gap of the bulk MgB 2 [jl6| . In tb 
second scenario the small gap is due to a weakened sr 
perconductivity on the very surface. Then two tunnelin 
contributions may exist: one corresponding to tunnelin 
to the weakly superconducting surface and another t 
the bulk MgB2 directly. Again, the feature at 7.5 m 1 
determines the energy of the bulk gap. 

A third and simpler interpretation is the intrinsic doi 
ble gap scenario for the bulk superconductivity in MgB; 
as predicted for the clean limit It originates from 
the existence of two different sheets of the Fermi sur- 
face, quasi-2D cylinders and a three-dimensional tubular 
network, and from a very particular anisotropic phonon 
spectrum. This results in a strong electron-phonon cou- 
pling with the quasi-2D bands and a much weaker cou- 



pling to the 3D band. Finally, two distinct and observ- 
able superconducting gaps appear in the quasiparticle 
spectrum, one larger and the other smaller than the BCS 
value. The effect is very sensitive to impurities, and van- 
ishes in the dirty limit in which just one gap exists. Tak- 
ing into account the very short coherence length £ < 50 
A |ll| , it is quite possible that both clean and dirty limit 
conditions may be found on the surface of such a granular 
sample. 

Further tunneling data are presented in Fig. 2 as strong 
evidence in favor of two gap superconductivity. Here the 
conductance spectra, selected from different locations of 
the same sample, clearly show double-peak structures of 
various relative amplitudes. Remarkably, the energy po- 
sitions of the gaps remain fixed, which is difficult to ex- 
plain in terms of a proximity effect. Indeed, in the latter 
case the DOS is a complicated function of the SN sand- 
wich parameters (interface transparency, metallic layer 
thickness, diffusion coefficient, etc.) and may be calcu- 
lated using Usadel's equations pCfl . In particular, both 
the apparent amplitudes of the peaks and their energy 
positions (especially those of the small induced gap) are 
very sensitive to any change of the SN sandwich parame- 
ters. This is evidently not the case in Fig. 2 (a), in which 
the relative amplitude of the gaps changes markedly, but 
their positions remain robust, clearly revealing two char- 

ftptpriQtir pnArcriPQ 




Sample bias, mV 



FIG. 2. (a) Tunneling conductance spectra (T = 4.2 K) 
normalized to unity and shifted for clarity (horizontal bars 
indicate zero conductance). The values of the gap widths 
(As and Al) are fixed but the relative strength of the peaks 
vary. The spectra can be fitted by only adjusting the coeffi- 
cient Cl the weight-factor of the large gap to the total density 
of states. The spectrum Cl = 0.13 was obtained in inverted 
SIN case 0. (b) Theoretical curves for the two-gap model, 
Eqs. (1-3), with different values for the tunneling probabilities 
Cl, Cs — ^-Cl- The case Cl=1 gives the usual BCS density 
of states with gap width Al, while Cl~0 (Cs = 1) represents 
the BCS case with gap width As- 



2 



The double-gap superconductivity explains naturally 
the above observations. While the peak positions are 
fixed, i.e. corresponding to two energy gaps, the spectral 
weight of each reflects the probability of tunneling to dif- 
ferent electronic bands (2D and 3D). Indeed, in powder 
based samples the orientation of grains with respect to 
the tunneling cone varies, as do the respective tunneling 
probabilities, from one location to another. This explains 
the variation of the relative amplitudes of the coherence 
peaks from grain to grain. As illustrated in Fig. 2 (b), 
the double-gap spectra can be generated by a weighted 
sum of two BCS density of states with coefficients Cl 
and Cs {Cl + Cs = 1) for the large and small gap DOS, 
respectively : 



a(V) = C L <J L (V) + C S v s \V) , 



(1) 



with 



r + co 

<j ls {V)k g(e-eV) p L < s (e)de, (2) 



/< s (e) = Re 



e-ir 



(3) 



Here g(e) = —df(e)/de is the usual thermal broadening 
function and T is Dyne's smoothing parameter. 

This simple analysis, Eqs. (1-3), appears sufficient to 
fit all of the double-gap conductance spectra of Fig. 2 
(a) (solid lines in the figure). We obtain statistically 
A L = 7.5 ± 0.5 meV and A s = 3.5 ± 0.4 meV with the 
smearing term not exceeding T — 0.2 meV. The double 
structure and other DOS features (zero bias conductance, 
spectral background, etc.) are perfectly reproduced by 
the theoretical function for different values of Cl- Even 
for the bottom spectrum in Fig. 2 (a), for which the exis- 
tence of the inner gap is not so evident, it was necessary 
to include both gap contributions. 

The double gap structure represents only a very small 
fraction of the acquired data over the surface. The STS 
technique allows to systematically explore different re- 
gions showing a variety of quasiparticle spectra. A com- 
plete I(V) spectrum is acquired locally at each pixel 
(x, y) of the STM image, so that a perfect correspon- 
dence between topographic and spectroscopic informa- 
tion is achieved. Moreover, a huge number of tunneling 
spectra are acquired (up to 65536 I(V) curves per image) 
allowing a significant statistics. The raw I(V) curves are 
numerically differentiated and presented either as a se- 
ries of tunneling conductance maps dl /dV(Vi, x,y), at 
a chosen bias Vi, or as tunneling conductance spectra 
dI/dV(V, Xi, yi) measured in particular locations (xi, yi). 
More details concerning our STS experimental setup can 
be found in S. 



In Fig. 3 we give an example of such an STS measure- 
ment. Image (a) is the topography acquired simultane- 
ously, of the surface (220 x 220 nm 2 ) showing relatively 
flat terraces, with a roughness of about 1-3 nm separated 
by steps 1-20 nm high. It does happen that the tip jumps 
from one region to another due to the surface roughness 
on a larger scale. A single occurence of this is indicated 
in (a) by the arrow. 
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FIG. 3. Complete DOS mapping of a 220nm x 220nm area 
of a chosen MgB2 grain, (a) Topographic image acquired in 
the constant current mode (sample bias V = -18 mV, I = 80 
pA). Flat regions are resolved. Black arrow indicates a tip 
jump during a scan, (b) Conductance map of the same area 
at V = -15 mV. (c) Conductance map at V = +5 mV. Bright 
and dark zones respectively indicate high and low conduc- 
tance zones, (d) In black: 'contaminated' regions (as found 
from (b) ); in white: selected regions (numbered from 1 to 8) 
described in the text. 

In Fig. 3 (b) we show the tunneling conductance map 
at V — —15 mV proving that the DOS is spatially in- 
homogeneous. Regions in which high-quality tunneling 
conditions are not met (the current is noisier, probably 
due to the surface contamination) are in black in Fig. 3 
(d). Of particular interest is the conductance map for 
V = +5 mV, i.e. at an energy close to the expected su- 
perconducting gap edge, Fig. 3 (c). In this map, bright 
and dark zones reveal the co-existing regions with high 
and low conductances. Comparing the topographic im- 
age (a) with the two conductance maps (b) and (c), we 
can easily select the most 'significant' regions in which 
noiseless tunneling conditions are achieved on a relatively 
flat surface. Such regions, of a typical size of 10-30 nm, 
are colored in white and numbered from 1 to 8 in Fig. 3 
(d). As a simple approach, the spectra are averaged over 
each selected region and the results are plotted in Fig. 4. 

The conductance spectra from the regions 1 and 2 are 
quite similar and both spectra present a smoothed re- 
duced gap of about 2.5 meV. The ratio 2A/k B T c = 1.48 
(Tc = 39 K for the bulk) is less than two times smaller 
than the weak coupling BCS ratio. Such spectra originate 
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either from the weakened Tc on the surface (due to the 
local surface degradation) or from a proximity induced 
gap. The conductance at zero-bias is far from zero and 
these spectra show a strong smearing term T = 0.3 A. 
Contrary to what was suggested in O], T does not re- 
flect the properties of the bulk MgB2 since it varies spa- 
tially, and even vanishes in other regions. This type of 
spectrum was also observed by Sharoni et al. Jl3| and 
explained in terms of the proximity effect. 
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ever the average gap value A = 4.6 meV, is still signif- 
icantly small, i.e. less than the ideal BCS value of 5.9 
meV by 20%. In view of the double gap spectra, as in 
Figs.l and 2, this characteristic single gap must corre- 
spond to the dirty limit. We conclude that about 99 % 
of the STS conductance spectra yield a single gap, with 
a probable reduced Tc, due to proximity effects or the 
dirty limit cases. Thus, the STS explains why the double 
gap DOS was not observed in previous tunneling reports. 

In conclusion, in this Letter we report the first direct 
evidence for two superconducting gaps in MgB2. These 
are due to different electron-phonon coupling in the two 
distinct parts of the Fermi surface, in agreement with 
the multiple gap model It is also consistent with 
the variation in relative spectral weight observed, due 
to the arbitrary orientation of the tunneling cone to the 
grains. The additional study of the temperature depen- 
dence of the two gaps supports our conclusions [[lTj . In 
the proximity effect or lowered surface Tc scenarios, the 
small gap should close at a lower temperature, while it 
should survive until 39 K together with the large gap 
in the double gap picture. No evidence for any large 
anisotropy of the order parameter nor pair breaking is 
found in the best spectra reported here. More signifi- 
cantly, the value of the large gap (7.5 ± 0.5 meV) leads 
to the ratio 2A/ksTc = 4.5 ±0.3 indicating the presence 
of strong electron-phonon coupling. STS performed on 
granular MgB2 at T = 4.2 K, shows the importance of 
spatial inhomogeneities in such samples. 

This work was supported by the Project ACI "Nanos- 
tructures" of the French Ministry of Research. 



FIG. 4. Normalized tunneling conductance spectra corre- 
sponding to the regions 1-8. The data points are shifted for 
clarity (horizontal bars indicate zero tunneling conductance). 
The spectra in the regions 1-2 present a small gap A = 2.5 
meV. These smoother spectra are fit (solid line) with a large 
r (0.3 A). In region 3, the gap width varies spatially, giving 
a smeared averaged characteristic, well fit (solid line) by two 
gaps of 4.3 mV and 5.7 mV and r< A. In the regions 4-6 the 
spectra have a BCS shape with a single gap. The theoretical 
fit (solid line) with A = 4.6 meV and F = meV is satis- 
factory. Regions 7-8: tunneling conductance shows no states, 
i.e. insulating behavior, in the energy window studied. 

The detailed analysis of region 3 (about 20 nm size) 
shows that two characteristic spectra are present: the 
first one, about 45% of the area, shows a gap of 4.3 
meV, while the second one has a gap of 5.7 meV (about 
55%). Thus, the average dl/dV curve in Fig. 4, having a 
small r = .05 meV, has a kink at the small gap en- 
ergy. Regions 4-6 are the most homogeneous and have 
the sharpest tunneling spectra with a clear single-gap 
BCS shape. The coherence peaks are very high and are 
situated at the same energy, and no states exist inside. 
Only thermal broadening is needed (T = meV), how- 
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